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\^ Abstract Using the analytical Fick-Jacobs approximation formalism 

,-H and extensive Brownian dynamics simulations we study particle trans- 

port through two-dimensional periodic channels with triangularly shaped 
walls. Directed motion is caused by the interplay of constant bias act- 
ing along the channel axis and a pressure-driven flow. In particular, 



we analyze the particle mobility and the effective diffusion coefficient. 

SThe mechanisms of entropic rectification is revealed in channels with a 
broken spatial reflection symmetry in presence of hydrodynamically en- 
-*— t forced entropic trapping. Due to the combined action of the forcing and 

the pressure-driven flow field, efficient rectification with a drastically 
CZ3 reduced diffusivity is achieved. 



1 Introduction 



Particle separation techniques of micro- or even nanosized particles are based on the 
fact that the particles' response to external stimuli, such as gradients or fields, de- 
pends on their size. Accordingly, conventional methods for filtering particles involve 
centrifugal fractionation pQ, phoretic forces [5] or external fields [5]. Recently, novel 
devices for particle separation were proposed |4|5|6)7j . These utilize the ratchet effect, 
i.e., the directed transport under non-equilibrium conditions in periodic systems with 
broken spatial symmetry |8l9j . The use of microfluidic channel systems is promising 
with respect to separation efficiency, speed and purity |5l7j . The transport in chan- 
nels with periodically varying walls exhibits peculiar transport phenomena |10lllll2] 
which can be treated by means of the so-termed Fick-Jacobs formalism and its gen- 
eralizations [13114115116117118119120] . The capability of such devices for separation of 
particles is rooted in the effect of entropic rectification, i.e., the rectification of motion 
^T caused by broken spatial symmetry |21|22|23|24j . 

Besides the direct forcing of the particle dynamics, the application of hydrody- 
namical flows presents an additional "degree of freedom" to control particle transport 
and optimize rectification. As pointed in a recent preliminary report by us [12] , upon 
combining a constant force causing the particle to move along the channel and a 
pressure-driven flow that drags the particle in the opposite direction results the phe- 
nomenon of hydrodynamically enforced entropic trapping was described; it implies 
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that for certain values of the constant force and the pressure drop, the particle mean 
flux vanishes and, in addition, the particles' diffusivity drastically reduces. With this 
work, we address a different problem, namely the topic of an entropic rectification in 
the presence of a hydrodynamic flow field. 

After presenting the model in Sect. [2] we apply the generalized Fick- Jacobs theory 
to the two-dimensional geometry with triangularly shaped confining walls in Sect.[3j 
We discuss the standard entropic rectification phenomenon in Sect. HI Sect. [5] is de- 
voted to the combined action of constant forcing and the hydrodynamical flow field 
on the rectification phenomenon. In Sect. [6] we summarize our main findings. 



2 Statement of the problem 

We consider spherical Brownian particles of radius R suspended in a solvent with 
dynamical viscosity r\ in a channel with symmetric confining walls of triangular shape, 
cf. Fig. IT] As argued in an earlier account |12j . we are interested in a two-dimensional 
(2D) geometry in which the planar channel is i-periodic, with the maximum and 
minimum widths Afi and Aw, respectively, and the position x s of the largest width. 
The confining zigzag walls are described by piecewise-linear boundary functions 

r\ ^ eL ilT + T L 6> ior x < x s (mod L) , 
UJ±{X)=± Y)^ + J ^, elsewhere. W 

^ X s — 1 1 — ' 

In Eq. (TTj) , 6 denotes the ratio of minimum to maximum channel width, i.e., 6 := 
Aw j AM, and e := (Afi — Alj) / L is the dimensionless geometric parameter |12ll9j . 
The channel's asymmetry is controlled by x s : for x s — > L/2 the channel exhibits 
reflection symmetry; for x s — > or x s — > L the largest asymmetry is achieved [24] , 
Note that for 8 — ¥ 1 and e —¥ the channel exhibits constant width Ail, w±(x) = 
±Af2/2. 

Assuming throughout (i) dilute particle density inside the channel, (ii) negligible 
particle size (R <C Aw), and (iii) a strong viscous dynamics, implies that inertial 
effects, hydrodynamic particle-particle and particle-wall interactions, and effects that 
can be initiated by rotation of particles can largely be neglected 12 25 26J. If so, the 
particles evolve in a laminar flow and the dynamics is well described by the over- 
damped Langevin equation. By measuring the lengths, energies, forces, and time in 
the scales of L, k^T, kgT/L, and the relaxation time r = 6irnRL 2 /(kgT), respec- 
tively, we obtain the dimensionless Langevin equation describing the particle's motion 

^=F(q) + V2€(t), F(q)=-Vtf(q) + Vx#(q), (2) 

where q = (x, y) T is the particle position, fc^T is the thermal energy, and the Gaussian 
random force | = (£x,£y) T obeys (&(<)) = 0, (&(f)£j(s) ) = S l:j 5{t-s); i,j G {x,y}. 
In general, any force field F exerted on particles can be decomposed into a curl- 
free part (scalar potential <P) and a divergence-free part (vector potential ?P r ), which 
constitute the two components of the Helmholtz's decomposition theorem. 

Hereafter, we study the impact of both contributions to the force F in Eq. (IT]) on 
the particle transport through the described geometry. On the one hand, we consider 
an external constant bias in x-direction with magnitude /, leading to ^(q) = —fx. 
On the other hand, we account for the difference between the particle velocity q and 
the local instantaneous velocity of the solvent u(q) based on the Stokes law, which 
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V(x*,y*)\ 
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Figure 1. A segment of a triangular shaped channel restricting the overdamped motion 
of point-like Brownian particles. The zigzag profiles with period L are given by periodic 
functions ui±(x), see Eq. dy, whose spatial asymmetry is determined by x s . The quantities 
Acu and AQ denote the minimal and maximal channel widths, respectively. The (blue) field 
lines depict a typical external force field F(q) exerted on the particles, which are the result 
of the interplay of a constant bias / and an oppositely oriented divergence-free force from 
the driven solvent induced by a pressure P(q) with pressure change AV along a unit cell. 
This interplay results both in vortices and stagnation points (red circles). 

gives us u(q) = V x tf'(q) with ^(q) = !f - (q)e z . Here, tf'(q) is the hydrodynamic 
stream function. As a result, Eq. (pi) turns into the Langevin equation 



dq 

df 



= fe x + u(q) + V2Z(t). 



(3) 



to be supplemented by no-flux boundary conditions for the particles at the walls. 



3 Generalized Fick-Jacobs approach 

We next present our generalized Fick-Jacobs approach PJ] which extends the standard 
Fick-Jacobs theory towards the most general force fields F as detailed with Eq. (pi). 
We first compute the joint probability density function (PDF) P (q, t) of finding the 
particle at the local position q at time t, provided it was at position q = (0, 0) at 
time t = 0. The evolution of P (q, t) is governed by the Smoluchowski equation |11I27| 



d t P(q, t) = -V [F(q)P(q, t)\ + V 2 P(q, t) 



(4) 



supplemented by the no-flux boundary conditions at the channel's walls. In the long 
time limit, lim^oo P (q, t) = P (q), the PDF has to satisfy the normalization condi- 
tion f . ,, P(q) d 2 q = 1, and be periodic, P(x + m,y) — P(x, y) , Vrn e Z. 

Assuming fast equilibration in the transverse channel direction, we perform an 
asymptotic perturbation analysis in the geometric parameter e -C 1 [19 20 , cf. Rcfs. 
[16117' 28], Upon re-scaling the transverse coordinate y — > ey, the profile functions 
and the vector potential become w±(x) — > eh±(x) and \[f —$ eH/, respectively. Ex- 
panding the joint PDF in a series in even powers of e, we get P(q, t) = Po(q, t) + 
e 2 Pi(q,t) + 0(e 4 ) and similarly <5(q) = # (q) + 0(e 2 ) and #(q) = %(q) + 0(e 2 ). 
Substituting this ansatz into Eq. Q and observing the boundary conditions, we ob- 
tain a hierarchic set of partial differential equations. For the steady state, we find 
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Po(q) = g{x) exp(— <Po(q)), where g(x) is obtained in the order 0(e 2 ). With the con- 
dition that the x-component of V x \P is periodic in x with unit period the stationary 

marginal PDF, P (x) = lim^oo f h 7^ dyP {q,t), yields 



Po(x) 



X- l I{x) 



(5) 



Here, I(a?) = e ^^ f* dx'e^^ >, with X = L dx I(x), and .F(x) is the generalized 
potential of mean force, reading 



F{x) = - In 



h + (x) 



dy e 



-*o(q) 



/i_(x) 



h+(x') 



|dx' /" dy (V x # ) a Pcq(yN') 



(6) 



h_(x') 



with P oq (y\x) = e *°( q ) / J, 'y'y dy e *o(<i). We reveal that F(x) comprises the usual 

entropic contribution (the logarithmic term) [10 29] caused by the non-holonomic 
constraint stemming from the boundaries [18 30 and an additional, energetic con- 
tribution stemming from &0 [T^j. In the absence of the flow field, i.e., for \P = 0, 
Eqs. (15]) and (pi) reduce to the commonly known result of the Fick- Jacobs approxi- 
mation [11114] . 

The kinetic equation for the time-dependent marginal PDF Po(x,t), with the 
steady-state solution in Eq. (p)]), is the generalized Fick- Jacobs equation 



Wt Po{x ' t) = dx- 



dF{x) 
dx 



P Q (x,t) 



J?! 

dx 



2 Po(x,t) 



(7) 



Based on Eq. ffi\, the stationary mean particle current for our problem reads [12] 



= km ^ = r 1 (1 - e Ar ) 



(8) 



wherein AT = J 7 (x+1)—J 7 (x). Note that Eq. pi) is a generalization of the well-known 
Stratonovich formula J31I which was originally derived for tilted periodic potentials. 



4 Entropic rectification 

Here, we shall focus on transport across triangular-shaped channel structures induced 
by a curl-free force field, P(q, t) = /, in a resting solvent, u(q, t) = 0. Accordingly, 
the generalized potential of mean force Eq. (fsj) simplifies to the "entropic" potential 
F(x)=A(x) = -fx-ln[2u(x)] \HM\- 

We first derive the dependence of the average particle velocity ( x ) on / and 
channel's parameters. According to Eq. ([8]), we obtain 



(1 



o /0^-m 



/ (l-x s )p(l-e-f) 



x s 6f 

e 1 ~ 5 



r 



x s Sf 



Jx B 



x s p{l-e-f) 



e 1 ~ s 



r[ 0, 



(x s 



1-5 



no, 



Xsf 



r[ 



■ 5 

1)/ 



(9) 



for the wall profiles Eq. (fTJ), where r (n, a) — f°° dte^ t t n ^ 1 is the upper incomplete 
gamma function. Another important transport quantity is the particle mobility \x 
reading 



ix :- 



l_x, 

f 



(10) 
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Figure 2. Numerical results (markers) for the particle mobility /i(/) (a) and mobility 
difference Afi(fo) = /x(/o) — /u(— /o) (b) versus the force magnitude / and /o, respectively, 
for a quiescent solvent, u = 0. For panel (a), the asymmetry parameter x 3 is varied, while 
the channel widths are fixed, Auj — 0.1 and AQ — 0.5. Note the logarithmic scale for the 
whole range of /. The horizontal dash-dotted line indicates unity. Inset (i): Comparison of 
Brownian dynamics simulations (markers) of Eq. Q with the Fick- Jacobs approximation 
(lines), Eq. Q9p, for a channel with AQ — 0.1 and Auj — 0.01. In panel (b), the aspect ratio 
S and the asymmetry parameter x 3 are varied, while AQ = 0.5. 



for any non-zero force /. In the case of infinitely large force strength /, the particle 
velocity converges to its free value, lim±f^ oc (x) — ±f, and therefore the particle 
mobility fj, tends to unity, regardless of the value of the asymmetry parameter x s G 
(0, 1). However, for a wedge-like shaped unit-cell, i.e., x s = or x s = 1, the particle 
mobility monotonically decreases with increasing |/| and converges to the asymptotic 
value /j,(f — » — oo ) = S for x s — and /i(/ — > oo) = S for x s = 1 [32] (not explicitly 
shown in Fig. [2]) . Referring to the Sutherland-Einstein relation, the effective diffusion 
coefficient -D G ff(/) '■— um t^oo(( x 2 ) — (x) )/2t coincides with the particle mobility 
in the diffusion dominated regime, i.e., /«1. For / — >• of Eq. ffih, we derive 



limMo(/) 



lim D eS (/) 
/->o 



2 (1 - S) 
(1 + 5) In (1/5) ' 



(11) 



Interestingly, within the presented Fick- Jacobs theory, the asymptotic value is only 
determined by the channel's aspect ratio i5 = Auj j ' AQ and, more importantly, it is 
independent of the asymmetry parameter x s . 

In Fig. [2] (a), we depict the dependences of the particle mobility /i on the mag- 
nitude of the bias / and its direction (positive and negative). The numerical results 
are obtained by Brownian dynamics simulation of Eq. pi) for 3 • 10 4 individual parti- 
cles. Note that there is a strong dependence on the direction of forcing 21 24 33 34J 
reflecting a broken symmetry of the channel (with respect to x). For x s 7^ 1/2, the 
channel's symmetry is broken, resulting in fi(f) ^ fi(—f). Starting out from a value 
smaller than the bulk value, i.e., fJ-(f) < 1, the nonlinear mobility monotonically in- 
creases with the strength of the bias |/| and approaches unity for |/| — > 00. However, 
for symmetric triangular channels corresponding to x s = 1/2, the particle mobility 
does not depend on the direction of the forcing and one observes /i(/) = /i(— /). 

If the force points into the direction of the gently rising flanks, i.e., for / < and 
x s > 0.5, the focusing character of the channel rules the propagation of the particle 
in channel direction while diffusion in transversal channel direction, which suppresses 
the mobility in channel direction, is reduced, cf. comparison of time scales in Ref. [29 . 
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Consequently, the nonlinear mobility monotonically increases with / and approaches 
the free value for / — > — oo, cf. Fig. [21(a). 

Interestingly, the dependence of /i on the forcing parameter exhibits a peculiar 
non-monotonic behavior if the forcing acts towards the direction of the steep flanks, 
/ > 0. Due to the steepness of the confining walls the force component directing 
the particles towards the bottleneck is rather small, ft oc (1 — x s ). Therefore the 
particles' residence time in one unit-cell increases, respectively, the particle mobility 
goes down with growing asymmetry parameter x s . Similar to the case of septate 
channels, x s — > 1, the particles can diffusively explore the full y-range which results 
in a pronounced minima of the nonlinear mobility for moderate forcing strengths. 
Nevertheless, as the focusing effect strengths more and more with increasing bias, the 
mobility tends to unity for / — > oo. 

In the inset of Fig. 0(a) we compare the numerical obtained results for /i (markers) 
with the Fick- Jacobs approach result (lines) , Eq. Q , for a weakly corrugated channel 
with AD = 0.1 and Alu = 0.01. Noteworthy, our exact analytic result matches very 
well the numerics for almost all values of / and x s . Deviations between Eq. ([9]) and 
the numerics occur only for extreme asymmetric channel geometries, i.e., for x a —¥ 1. 
In contrast to the usually studied sinusoidally modulated channel [12I19I29I35J . the 
requirement £ <§; 1 is not sufficient to ensure the validity of the Fick- Jacobs approxi- 
mation (leading order in e). Concurrently, as an analysis of the different involved time 
scales predicts, the cross-section variation rate A = (1 — x s )/e [SS] has to be small 



as well. Note that in contradiction to Eq. (II) (zeroth order perturbation theory), 
our numerical simulations show a weak dependence of lim/_>o /•*(/) on the asymmetry 
parameter x s . 

We next discuss the impact of the channel's asymmetry, 0.5 < x s < 1, on 
the entropic rectification properties. We assume that the external bias instantly 
switches between the two values ±/o after half a period T, i.e., /(£) = ±/o. Sim- 
ilar to the Brownian motor phenomenon [489, such oscillating forcing results in 
directed motion along the a;-axis. Assuming adiabatic driving, i.e., the period is much 
longer than any other involved time scale, the effective drift velocity is given by 
((x)) T :— 1/T L dt (£(/(£))) and reduces for dichotomic, unbiased driving to 

(Wo)»r= ^((i(/o)) + (i(-/o)))= ^m(/o)/o, (12) 

where Afi(fo) := /u(/o) — m(~ ' /o) denotes the mobility difference for forward and back- 
ward driving. Note that the relative mobility difference a = |Z\/i(/o)|/ (m(/o) + m(~ /o)) 
for quantifying the rectification was introduced earlier as another measure for entropic 
rectification [2Tj . 

In Fig. 0(b), we present Z\/x(/o) for different channel asymmetries x s and aspect 
ratios 6. While for symmetric channels, x s = 0.5, the mobility difference equals zero 
regardless the value of /o, A[i(f ) — 0, a non-monotonic behavior is observed if the 
channel symmetry is broken. For 0.5 < x s < I, as fx(-fo) > A*(/o) the mobility 
difference A[i is negative and thus the particles' motion is rectified, viz. they move 
with ( ( x) ) T < to the left. For the limiting cases f — > oo and f — > 0, the forward 
and backward mobility equals each other, cf. Fig. [2] (a), and correspondingly A^(f ) 
approaches zero. Note that both the minimum value of Afj,(fo) decreases and the 
position of the minimum shifts to larger force magnitudes /q with growing channel 
asymmetry. Solely for wedge-like shaped unit-cells, x s = 1, Afj,(fo) monotonically 
decreases with /o from zero to the asymptotic value Afi(oo) = (6 — I). Decreasing 
the bottleneck width, respectively, the aspect ratio 5 results in a systematic shift of 
all curves towards smaller values for A/j,(f ). Simultaneously, the value of /o where 
the mobility difference is maximal remains. 
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As different sized particles may be treated by consideration of an effective bottle- 
neck widths |37j and given that the mobility difference exhibits a strong dependence 
on 5, our channel setup is an ideal candidate for a device separating particles according 
to their sizes [5]. Since the particle transport is accompanied by a strong enhance- 
ment of the effective diffusivity [35] , the quality of entropic particle rectification may 
be weak. Recently, we found that the counteraction of a pressure-driven flow and a 
constant bias of strength / ensues the intriguing hydrodynamically enforced entropic 
trapping (HEET) phenomenon where the vanishing of the mean particle current is 
accompanied by a significant suppression of diffusion [12] . Consequently, the quality 
of hydrodynamically induced particle rectification may be better. 



5 Interplay of pressure-driven solvent flow and external bias 

Next, we reveal the role of a flow field u(q) = V x tf'(q) on the particle transport 
through the triangular channel structures, presenting an example of the divergence- 
free force. The dynamics of particles is described by Eq. ^, which implies an one-way 
coupling between the solvent and the particles; for a dilute suspension, the fluid flow 
influences the particle dynamics but not vice versa [38j . 

A slow viscous steady flow of an incompressible solvent is governed by the dimen- 
sionless Stokes or "creeping flow" equations |25l39j . 

VP(q) = V 2 u(q), V-u(q) = 0, (13) 

being valid for small Reynolds number Re = pL 2 /(r]T) <C 1. The flow velocity 
u= (u x ,u y ) and the pressure 'P(q) are measured in the units of L/t and tj/t, 
respectively. Further, we require that u obeys periodicity, \i{x.y) — u(x + l,y), and 
the no-slip boundary conditions, u(q) = 0, Vq € channel wall. The pressure satisfies 
V{x + 1, y) = V(x, y) + AT, where AV is the pressure drop along one unit cell. 



Applying the curl to both sides of first relation in Eq. ( 13 ) eliminates 'P(q), yielding 
the biharmonic equation V &(x, y) = for the stream function &, related to the flow 
velocities as u x = d y ty and u y — ~d x W. Following the asymptotic procedure described 
in Sect, pi we solve the biharmonic equation = dt&o(x, y) + 0(e 2 ) with the no-slip 
boundary conditions, d y Wo — at the channel walls y = h±(x), and the conditions 
specifying the flow throughput Q = -AV/(12(n~ 3 {x)) x ) [3D], W = at y = h-(x) 
and $0 = Q at y = h + (x), to obtain in leading order in e: 

AV \y- h_(x)f [3h + (x) - h-(x) - 2y] 
12 H»(x) {H-»(x)) a ' l ' 

where H(x) — h+(x) — h-(x) is the re-scaled local width and (-) x — J Q • dx denotes 
the average over one period of the channel. 



5.1 Purely flow-driven transport 

In the absence of conservative forces, when <?(q) = (/ = 0), and for no-slip boundary 
conditions at the walls, the joint probability density function is spatially uniform, 
P(q) = 1/ {%{x)) , unless hydrodynamic interactions between particles and walls 
come into play |25|41j . 

Then, for the purely flow-driven case, / = 0, the mean particle current (i) ° w = 
CH{x))~ f Q dx J h + M dyv,Q(x,y) is explicitly evaluated for the wall profiles given by 
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Figure 3. The scaled mean particle velocity ( x } (a) and the effective diffusivity D c g (b) as 
functions of the pressure drop AV for a constant bias / = 100 and different asymmetry pa- 
rameters x„ for triangularly shaped channels with maximum and minimum widths Aft = 0.5 
and Aoj = 0.1, respectively. The dash-dotted line in panel (a) indicates the analytic estimate 
for / = 0, Eq. (15 I. The horizontal line in panel (b) corresponds to the value of the bulk 
diffusivity. The markers are for Brownian dynamics simulations of Eq. pi with the flow field 
obtained via numerically solving Eq. ( |13[ ). Inset (i): The dependence of critical force value 
/or on the critical pressure drop AVer for different aspect ratios S is depicted, viz. 8 — 0.1 
(blank markers) and S = 0.2 (filled markers). The analytic result Eq. (16 1 is shown by the 
dashed lines. 



Eq. |l]) to give 



\ flow 



AV 



AP{ADf 2<5 2 ln(l/<5) 



i2<«>„ n- 



12 



1-S 2 



(/ = 0). 



(15) 



The mean particle current is as for the Poiseuille flow between plane parallel walls, 
(x) s=1 = -AV{AQ) 2 /12, modified by the factor k{8) := 2S 2 In {1/5) / (l - S 2 ) ac- 
counting for the corrugation; note that lim^i k(6) = 1. Thus, the solvent flows from 
left to right for AV < and vice versa for AV > 0. Remarkably, result Eq. ( [15] ) 
is independent of the asymmetry parameter x s . Moreover, since (x(AV)) Q = 
— (x(—AV) ) ° w , no rectification of point-sized objects can be achieved by periodi- 
cally switching of the flow's direction - even for asymmetric channels. This result, 
however, has no contradiction with the "drift ratchet" [7], where the finite size of the 
particles was crucial for the transport. 



5.2 Interplay of solvent flow and external bias 



In Fig. Km a), we depict the mean particle velocity (x) as a function of the pressure 
drop AV for a constant bias / and different values of the asy mmetry parameter x s . 
While ( x ) exhibits the linear dependence predicted by Eq. ( 15 ) for / = 0, its behavior 
changes drastically for / ^ 0. For \AV\ 3> 1 the drag force dominates over the external 
bias and thus (x) oc —AV. For intermediate values of AV, (x) is nearly constant 
in a broad range of \AV\. For AV > 0, the constant bias drives the particles in the 
direction opposite to the flow. At a critical pressure drop AV cr , a sharp transition 
from positive to negative values of (x) is observed, cf. Ref. [12J . As follows from 
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Eq. Hi, (£) =0 when F(x + 1) 
\AV 



1 
12 



Fix) 
(U(x 



AT 



[H(x) 



-3' 



0, yielding for the critical ratio 

(Auj) 2 ln(l/5) 
' 6(1 - J 2 ) ' 



(16) 



being independent of asymmetry parameter x s , cf. Fig. p5[i) . 

The described behavior of (x) is reflected by peculiar features in the effective 
diffusion coefficient, cf. Fig. [3lb). If for \AV\ <C 1, -D e ff is mainly determined by the 
channel's geometry and the constant bias /. At AV sa AVer, a drastic reduction of 
the diffusivity can be observed. This effect, which occurs when the constant bias and 
the flow start to counteract such that the field F(q) = fe x + u(q) contains vortices 
and stagnation points and (x) ~ 0, is referred to as the hydrodynamically enforced 
entropic trapping (HEET). With the further growth in AV, D c g exhibits Taylor- Aris 
dispersion [42I43J irrespective of the channel constriction, i.e., D g oc (Af2 (x)) /192. 
Thus, the interaction with the flow can significantly affect the macroscopic quantity 
D cS . 
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Figure 4. Scaled effective drift velocity {{x)) T (a) and effective Peclet number 
{{x)) T / {D a g) T (b) for various compartment shapes. The Brownian particles are driven 
by a constant pressure drop AV = 10 4 and a dichotomic switching bias f(t) = ±/o. The 
horizontal dashed line indicates unity, whereas the vertical dash-dotted lines show the value 
of / cr , cf. Eq. ( 16 1. The channel parameter are set to AQ — 0.5 and Auj — 0.1. 



The strong variation of the mean particle velocity ( x ) accompanied by the drastic 
reduction of diffusivity D c g at the critical conditions enables the possibility for effi- 
cient rectification. In Fig. HI we present the impact of the asymmetry parameter x s on 
the effective drift velocity ( ( x ) ) T and the effective Peclet number ( ( x) ) T / ( D a g ) T 
for a dichotomic switching bias f(t) = ±/o. Since the mean particle velocity is domi- 
nated by the Stokes' drag force for / < / cr the time-averaged (over period T) velocity 
( (x) ) T coincides with the undisturbed particle velocity (x(AV) ). As f approaches 
/cr, ((%))t g rows an d thus ((x)) T > (x(AV)) independent of x s . Interestingly, 
((x)) T > (x(AV)) remains valid for strong bias /o despite that the mean value 
of the dichotomic force equals zero, (f(i)) T = 0. Thereby the enhancement of (x) 
grows with x s . We emphasize that enhancement of the particle transport is accom- 
panied by a resonance-like behavior of the effective Peclet number ( ( x) ) T / ( D CT j ) T , 
see Fig. Wlb). We find that the latter attains a local maximum at the critical force 
magnitude / = /„■, whose value is given by (x(-f ,AV) ) /D eS (-f ,AV). While 
the effective Peclet number grows with / for a symmetric compartment, x s — 0.5, it 
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goes to zero for an asymmetric one, see x s = 0.99, as caused by the strong enhance- 
ment of the effective diffusivity in triangularly shaped channels [36 . For / — > oo, we 
expect that ( ( x ) ) T / ( D c g ) T starts to grow since D e ff converges to unity for infinite 
strong external bias. 



6 Conclusions 

In summary, we addressed the problem of transport of Brownian particles in a two- 
dimensional generally asymmetric channel with periodic confining walls of triangular, 
zigzag profile. To study rectification mechanisms, we investigated the interplay be- 
tween a constant forcing acting on the particles along the channel axis and the Stokes' 
drag stemming from a pressure-driven flow field. Along with pure entropic rectifica- 
tion observed for channels with a broken spatial reflection symmetry in the presence 
of a constant external bias, the hydrodynamic flow field paves the way to efficiently 
rectify the particle transport utilizing the effect of hydrodynamically enforced entropic 
trapping, being accompanied by a drastically reduced diffusivity. While the usual en- 
tropic rectification phenomenon was proposed as the basic mechanism for efficient 
separation of particles by generating opposite fluxes for particles of different sizes [5] , 
hydrodynamic flows can be used to achieve lower diffusivity of particles resulting in 
an increased narrowing of distributions for the target particles. The latter feature 
may prove advantageous in tailoring more efficient separation devices as compared to 
currently used methods in separating micro- and/or nanosized particles. 
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